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The article below reproduces, with some minor changes, the text of the inaugural lecture 
delivered by the author on 6th February 1961 as professor extraordinary of botany at the Univer- 
sity of Utrecht. The author sketches the “balance of nature’’ as a state of equilibrium in con- 
stant evolution, from primeval forest to single-crop farming, and from the attendant menace of 
agricultural pests to the present practice of combating these pests with chemical agents. He 
draws attention to the curious fact that plants themselves use chemical agents in an interne- 
cine struggle for existence, and develops the intriguing concept of a “negative biology’’ based 
on the action of pesticides, which promises to penetrate far into the secrets of living matter. 

The text published here is supplemented, with the kind assistance of the author, by illustra- 
tions, a bibliography and an appendix of structural formulae. 


Many publications addressed to the general public 
speak of the “balance of nature”. They may often 
give the impression that this balance is more or less 
labile, so easily does it seem liable te disturbance. 
One might imagine that nature is comparable to 
a chemical balance, so that a slight shift in the 
balancing forces implied might be expected to 
“tip the scales” right over. 

An equilibrium as achieved on a chemical balance 
is however never encountered in nature. Darwin ') 
and others have described at considerable length 
how the balance of nature in fact operates. 


The primeval forest 


Darwin was one of the first to see the tropical 
forest not as a community of plants in peaceful co- 
existence, but on the contrary as the scene of a 
silent but merciless struggle for food, space and, 
above all, light. The weapons employed in that 
struggle are extraordinarily varied. Some plants 
exploit their physiological aptitude for vigorous 
upward growth to break out rapidly above the gloom 
of the lower levels and so offer a broad crown to the 
light. Other plants are so physiologically adapted 
as to be able to populate these lower levels ( fig. 1). 

I mention the word physiological because, as a 


physiologist myself, I think in the first place of such 
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properties as the highly economical photosynthesis 
required for the formation of enough assimilates at 
low light intensities; of growth hormones, like 
auxins and gibberellins, which induce exceptional 
elongation in many plants; and of factors which, 
in spite of that elongation, enable the plant to 
develop a sufficiently sturdy stem. 

There are other plants, e.g. lianas, whose excessive 
elongation prevents them from acquiring the neces- 
sary sturdiness, and yet they still contrive to reach 
the higher levels by climbing up other planis, which 
may ultimately collapse under their burden. 

Others again begin their lives in the light by ger- 
minating on an upper branch of a fully-grown tree. 
Such plants are called epiphytes. An epiphyte has 
to be extremely sparing with the mineral salts avail- 
able to it, for it cannot in general draw them from 
the soil. If it is capable of sending out aerial roots, 
as some species of Ficus do, the roots develop down- 
wards and may ultimately reach the ground and 
penetrate its surface; a tree is thus produced which, 
as it were, grows from top to bottom. The aerial 
roots multiply and become rapidly thicker. The epi- 
phyte gradually proliferates around its“ host”, and 
at last destroys it. The “tree strangler’’, now firmly 
anchored by the lignification of its original aerial 
roots, spreads out its crown in the place once occu- 
pied by its host. 

Although the innumerable seeds that fall to the 
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floor of the forest may well germinate, their chance 
of survival is remote. The struggle is ceaseless and 
unrelenting, and only the best adapted and most 
fortunate can develop to maturity. But I need 
dwell no further on this fierce competition for 
survival in the primeval forest; it is a subject with 
which we have grown sufficiently familiar in the 
hundred years and more that have elapsed since 
the appearance of Darwin’s brilliant work ““On the 
origin of species”. 

The same fierce struggle as in the primeval forest 
is fought in practically all plant and animal com- 
munities. Here too the balance of nature is quite 


unlike a chemical succeeds in developing to maturity 
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of the difficulties which are involved in agricul- 
ture, in eradicating the existing plant growth 
and cultivating useful species. The surrounding 
vegetation will have constantly fought back to 
recover its lost ground. Only by dint of hard work, 
ceaseless weeding and the most intensive cultivation 
of the soil can they ever have succeeded in raising 
their plants to maturity. Not for nothing is it written 
in the old Dutch “Katechismus der Natuur” that 
God created weeds besides the useful herbs, for 
had He not done so, the farmers would have grown 
idle *). 

Plainly, then, the agricultural population is inter- 


fering with the balance of nature in a singularly 


Fig. 1. Vegetation in the primeval forest: Ruwenzori massif in the Congo. The photograph 
was taken in a fairly open spot. The characteristic struggle in which the plants deprive 
each other of essential sunlight is naturally difficult to represent in a photograph. 


disturbs the existing pattern of relationships to some 
extent, for it intrudes itself among the other or- 
ganisms, often destroying its neighbours, and in so 
doing changes the order of things. By the very fact 
of preserving itself, every living being, whether 
plant or animal, has a considerable influence on its 
environment, which would certainly have looked 
different without it. Thus, it can be said that the 
balance of nature is constantly being disturbed to 


some extent by each plant and each animal. 


Man and the balance of nature 


Man too is a part of nature in this context, and 
to preserve himself is bound to shift the established 
balance, otherwise there would be no room for him. 
The early tillers of the soil must have been only 
balance. Every that too 


organism well aware 


cruel fashion in felling trees and burning woodland, 
in ploughing and otherwise mechanically working 
the land. With the growth of world population this 
interference is assuming enormous proportions. In 
more densely populated regions little remains of the 
original flora and fauna. Even so, the best adapted 
species have lost nothing of their aggressiveness, 
and the farmer still has no chance to be idle. 

The practice of agriculture requires that large 
acreages be devoted to the specialized cultivation of 
individual crops (single-crop farming). As a logical 
result, adapted species amongst the other organisms, 
whether weeds, insects, eelworms or moulds, have 
been able to make large-scale invasions (figs 2and 3). 

This brings us to the diseases and pests of agricul- 
ture. Since almost every animal and plant has its 


own peculiar parasites, and every agricultural plague 
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Fig. 2. Land reclaimed from inland seas may often be invaded by adapted weeds. For 
example, during the drainage of the 140 000-acre Oost-Flevoland polder, formerly a 
part of the Zuijder Zee, broad expanses were conquered by “marsh endive’’ (Senecio pa- 


luster). 


is caused by an animal or plant, pests which flourish 
in a given cultivated area are always followed by 
their own parasites. 

The demands made on agriculture are steadily 
growing, and the damage caused by diseases and 
pests are felt more and more keenly. It is therefore 
not surprising that the agriculturist must take 


increasingly drastic measures to derive full benefit 
from his labour. 

People have been so thoroughly accustomed to 
modern methods of farming, and to the periodic 
recurrence of pests and diseases, that they regard this 
state of things as a natural balance. If, for example, 
after having exterminated large numbers of insects 


Fig. 3. Another part of the polder shown in fig. 2, in a more advanced stage. The weeds 
which were formerly present have been replaced by rapeseed. (The 86000 acres under 
n this polder form the largest area of its kind in Europe.) In the meantime, 
another pest has developed in this single-crop area: the fields of rapeseed have become 


thickly infested with weevils. 


rapeseed i 
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in his orchards by spraying with DDT [1] *), the 
fruit-grower finds that his action has enabled the 
mite population to multiply, he promptly blames 
himself for having disturbed the balance of nature. 
Nothing is farther from the truth; that balance was 
disturbed long ago, and now it has only been shifted 
a little further. 

Usually it is the chemical control agents that are 
held responsible for such disturbances of the balance 
of nature. The physical methods, such as felling and 
burning, tilling, drainage, irrigation and the like, 
are generally accepted without demur. 

It need hardly be said that there are no real 
grounds whatever for making this distinction be- 
tween “inhumane chemical’ and “humane physi- 
cal’’ methods; it springs simply from the mistrust 
of innovations and the concomitant idealization of 
the old and familiar. 

In fact, nature itself makes widespread use of 
chemical means of defence and aggression. Many 
plants, for instance, protect themselves from insects 
by toxic alkaloids. I need only mention the upas tree 
of South-East Asia (Antiaris toxicaria), whose milky 
sap is used as arrow poison. 

Another mode of chemical warfare, this time be- 
tween plants among themselves, takes place below 
ground. The roots of some plants excrete substances 
which, if absorbed by the roots of other plants, 
severely damage the latter or inhibit their growth. 

This form of root competition is much more com- 
mon than might be supposed. Many think of root 
competition as a struggle for mineral salts or water. 
That is by no means always the case. Bonner dis- 
covered that the rubber-producing guayule plant 
secretes cinnamic acid [2] from its roots, which acts 
as a growth-hormone antagonist when taken up by 
other plants and so inhibits their development *). 
In Java we were able to isolate from the soil in which 
Salvia occidentalis was growing a substance capable 
of almost completely inhibiting the growth of young 
coffee plants *) ( fig. 4). 

It was reported some years ago that when land 
badly infested by eelworms is planted with African 
marigolds (Tagetes), the eelworm population rapidly 
declines. This observation attracted considerabbe 
attention because of the difficulty of controlling 
this pest effectively. Investigations by Uhlenbroek 
and Bijlo showed that the roots of Tagetes contain 
a substance, terthienyl [3], which is extremely 
toxic to eelworms °). These are just a few of the 
many indications that the underground struggle in 


nature is fought, at least in part, by chemical means. 


*) The figures between square brackets refer to the structural 
formulae in the appendix. 
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Fig. 4. Example of root competition. The coffee plant on the 
left was grown together with Salvia occidentalis in a water cul- 
ture (after which the Salvia was removed). The coffee plant on 
the right was grown without Salvia. The growth of the left- 
hand specimen has been severely inhibited by substances 
excreted from the Salvia roots. 


Above the ground, too, nature makes frequent 
use of chemical agents. One example is the shrub 
Encelia which has been studied by Went °) land 
Bonner *). This plant, a native of the Californian 
desert, exudes from its leaves a chemical substance 
that inhibits the germination of seeds. The sporadic 
rainfall in these urid regions washes the substance 
on to the ground under the shrubs. The seeds lying 
there would normally germinate in these humid 
intervals, but the inhibitor washed from the leaves 
prevents them from doing so. In this way the En- 


celia wards off the competition of other plants that 
might otherwise overwhelm it. 


Chemical control of pests, a necessary measure 


To control the spontaneous spread of weeds, the 
chemical industry has developed numerous _her- 
bicides which, administered in very small quantities, 
apparently block an essential biochemical reaction 
in plants. 

Natural vegetation, already so impoverished by 
mechanical methods of agriculture, is menaced even 
more by these chemical agents. The botanist, of 


———————— 
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course, views this development with regret. It is 
appalling to imagine what the Alpine meadows would 
look like if weed control of this kind were rigorously 
applied there. Our own native pastures have lost 
much of their former beauty now that they are no 
longer bright with buttercups, ladysmock, dan- 
delions and sorrel, not to mention the wild orchid 
and other plants now rarely to be found. 

We feel this as a sad deprivation and forget that 
the severest impoverishment was caused long ago by 
the growth of the population, resulting in the de- 
struction of forests and the introduction of mechani- 
cal methods of cultivation. Chemical methods are 
simply a further step forward in a development that 
was initiated by our distant forebears, a develop- 
ment which is making the creation of nature reserves 
more and more urgently necessary. 

Our most sensible policy is undoubtedly to take 
a realistic view of these chemical innovations in 
agriculture, and to set bounds to their application 
by establishing reserves where their use is forbidden. 
This whole course of events when all is said and done 
is the consequence of the enormous expansion of 
world population, which is a much more disquieting 
phenomenon. 


Biochemical equilibria in cells 


Before dealing with these pesticides any further, 
I should like to look for a moment at the structure 
of living matter, at the organism itself. There is a 
natural balance in an organism too, but of quite a 
different kind. Here it is a complex interplay of 
biochemical processes in which the chemical re- 
actions are on the whole reversible. Each reaction 
thus represents a state of dynamic equilibrium. If 
the product of a reaction is removed by a subsequent 
reaction, the first reaction will continue to proceed 
in the same direction. If the reaction product accu- 
mulates, however, the reaction will come to a stand- 
still or may even be reversed. 

Nearly every compound formed in an organism 
should be seen as a link in a long chain which, with 
the uptake of energy, can lead to more intricate 
compounds, or which, upon the dissipation of en- 
ergy, is normally degraded to simpler compounds. 

The picture of a compound as a link in a long 
chain is really oversimplified. In fact we should 
think of it as a node in a network, indeed in a three- 
dimensional network. A system built up in this way 
can clearly be a stable one. Again, many vital pro- 
cesses are stabilized by ingenious feedback mech- 
anisms. If a line is disturbed anywhere, the process 
usually continues by a roundabout route. Only the 
blockage of the most essential lines results in the 
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death of the organism. Now what are these essential 
lines? 

Indications in this direction may well be provided 
by the multifarious herbicides, fungicides and in- 
secticides in use, for it is precisely these substances, 
applied in such extremely small concentrations, that 
yield such remarkable results, albeit in terms of the 
destruction of living matter. The study of these 
agents thus constitutes a sort of “negative biology”. 
The following example will serve to illustrate the 
value of this negative biology in the fields of physi- 
ology and biochemistry. 

A line of major importance is that along which 
energy is supplied in living matter. The energy 
required for synthesis is derived in most cases from 
the compound adenosine triphosphate [4], generally 
abbreviated to ATP. This compound is formed in 
processes of photosynthesis (photosynthetic phos- 
phorylation), respiration (oxidative phosphory- 
lation) and fermentation “). 

In all three processes the energy made available 
by the oxidation of reduced compounds is stored in 
a phosphate bond of the ATP, which can then supply 
the energy for further processes that need it. 

Any substance that impedes the formation of ATP 
must therefore be highly toxic to the organism. In 
agriculture substances of this kind are used which 
inhibit ATP formation in the process of oxidative 
phosphorylation. Known as “uncouplers’”’, these 
agents do not stop respiration — on the contrary, 
they speed it up considerably — but no ATP is 
produced. The dinitrophenols are powerful uncoup- 
lers, and are used in large quantities as herbicides 
and insecticides (DNBP, DNOC [5}). 

The dinitrophenols are also used as respiration 
uncouplers in many scientific investigations. Proces- 
ses requiring energy are stopped by such a sub- 
stance, as the necessary energy is normally supplied 
to the process via the ATP. In this way it is possible 
to distinguish energy-consuming processes from 
those that proceed freely. The dinitrophenols have 
been used for many years in physiological and 
biochemical research on e.g. the uptake and trans- 
location of substances in plants. A good deal has 
thus come to be known about their action. 

Another weed killer well known to agriculturists, 
2,4-dichlorophenoxyacetic acid [6], belongs to the 
group of growth substances, its action on the plant 
being analogous to that of the normal growth hor- 
mone, indoleacetic acid [7]. Whereas, however, indo- 
leacetic acid is a link in a chain of enzymatic re- 
actions and is therefore dependent on other reactions 
both for its synthesis and for its breakdown, the 
dichlorophenoxyacetic acids do not fit into the 
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biochemical scheme of the plant. They are assimi- 
lated and function as if they were growth hormones, 
but they cannot be broken down in the normal cycle 
and as a result the growth and development of the 
plant are completely disorganized ( fig. 5). 

There are some other chemical agents I should 
now like to mention because of their potential value 
to fundamental scientific research. Indeed, through 
their apparent ability to disturb the biological bal- 
ance inside the cell, they may well provide us with 
the key to the essential processes in living matter. 
Unfortunately I must add at once that in this respect 
they have received far too little attention. What is 
at present known of their physiological and biochem- 
ical action stems mainly from the laboratories of 
the industries that produce them. 

Years after the introduction of chloropheny]- 
dimethyl urea [8], called “Monuron”, as a herbicide 
in agriculture, nothing was known about the mech- 
anism of its action. Wessels concluded on theoretical 
grounds that it was probably a potent inhibitor of 
photosynthesis 5). This in fact proved to be the case, 
further research showing it to be one of the most 
active of such substances. So specific an inhibitor 
will undoubtedly be of considerable value to research 
into the process of photosynthesis. 

Then there are the substituted triazines which, 
under the names of “Simazin” [9] and “Atrazin”, 
are used in large quantities as weed-killers. It is 
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known that these also inhibit photosynthesis, but 
it is quite likely that they arrest other vital processes 
too. 

These are just some examples of substances about 
whose active mechanism we have at least some indica- 
tions. But there are numerous other herbicides whose 
active mechanism is still completely unknown. One 
such substance is the germination inhibitor 2,6-di- 
chlorobenzonitrile [10], minute doses of which make 
it impossible for most seeds to germinate, and more- 
over completely stop the growth of young plants 
( fig. 6). Others in this category are marketed under 
the proprietary names “IPC”’,“Avadex”’, “Dalapon”’, 
“Amitrol” and “Karsil”. Further, there is “Reg- 
Jone”, a bromodipyridyl compound, extremely small 
concentrations of which (less than 1 pound per acre) 
are sufficient, in some way still unknown, to derange 
the physiological balance in plants so thoroughly 
that they die off in a very short time. It would be 
highly interesting to discover the particular system 
on which such minute quantities act. 

Amongst the many fungicides in agricultural use, 
the only ones whose active mechanism is not com- 
pletely unknown are the dithiocarbaminates [11]. 
Their physiological action has been elucidated, at 
least in part, by Kaars Sijpesteyn and co-workers °), 
who also found the explanation for the double peaks 
in the curves showing the relation between effect and 


concentration ( fig. 7). 


f ig. 5. The effect of the herbicide 2,4,5-trichlorophenoxyacetic acid on brambles is seen 
from the treated shrub on the right. (The shrub on the left is untreated.) This herbicid 
acts as a growth hormone, but does not fit into the biochemical scheme of the pl : 
and therefore completely disorganizes its development. j acs 
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Of the insecticides the group of phosphorus com- 
pounds is the subject of a great deal of biochemical 
research, owing to the fact that they inhibit the 
enzyme choline-esterase. Since the latter is essential 
to the proper functioning of the nervous system, the 
medical profession is keenly interested and various 
medicophysiological and biochemical laboratories 
are at present studying these substances. 

Apart from this group, there are many more in- 
secticides and acaricides which have so far only been 


studied as regards their practical application. 


Fig. 7. Example of complications in the action of chemical 
control agents. The curves show the relation between the con- 
centration of NaDDC [11] and its effect on Aspergillus niger 
in the presence and absence of copper (white and black circles, 
respectively). Where copper is present, the curves show two 
peaks. An explanation of this phenomenon will be found in 
reference °). The ordinate gives the content p of accumulated 
pyruvic acid, which is a measure of the biochemical disturbance 
on which the action of this pesticide depends. 


Fig. 6. The effect of 2,6-dichlo- 
robenzonitrile (code number 
H{ 133), in various concentra- 
tions, on a) germinating millet 
(agar-agar culture) and b) young 
potted oats. In b) the agent 
was applied in vapour form 
(left and centre in  concen- 
trations of 0.3 and 0.1 mg/litre 
respectively, right untreated). 
In c) the effect of the agent asa 
weedkiller was tested on a field 
of gladioli. The part in the 
foreground was sprayed with 


1.6 kg/acre of H 133. 


Conclusion 
It is my opinion that there are many chemical 


agents used in agriculture that are certainly worth 


intensive and fundamental research to discover 
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the mechanism of their action. Some of these sub- 
stances are already subjects of such research. It is to 
be hoped that many of the others will soon be sub- 
jected to a similar investigation. The best way to 
make research along these lines yield rich fruit 


apart from properly equipping the laboratories in 
question — is close cooperation between biologists 
on the one hand and biochemists and agriculturists 


on the other. 


Appendix: structural formulae of some compounds mentioned 
in the text 
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Summary. Principal contents of the inaugural lecture delivered 
by the author as professor extraordinary of botany at the 
University of Utrecht. In e.g. cultivating the land, man par- 
ticipates in the constant evolution of the “balance of nature’. 
The chemical agents used to this end should not be seen merely 
as a means of pest control; they can also be of value to the 
fundamental study of living matter. Their action relies on the 
disturbance of essential processes in the living organism, and 
because of that fact they may well be the key to understanding 
those processes. The author advocates close cooperation between 
biologists, biochemists and agriculturists with a view to 
more intensive research into this “negative biology”, and 


mentions various chemical agents that might be studied along 
these lines, 
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THE APPLICATION OF CONTROL THEORY TO LINEAR CONTROL SYSTEMS 


by M. van TOL m)2 


621-53.001 


Up to the beginning of the ’forties, control engineering was approached largely on an em- 
pirical basis; measures taken to improve a control system tended to rely more on practical 
experience and intuition than on theoretical insight. Improvement came when it was recognized 
that the theory of control systems and servomechanisms could be derived directly from the theory 


of negative-feedback amplifiers. 


The article below discusses the present-day treatment of control problems with the aid of 
complex variables, a method commonly used in alternating-current theory. As one example, 


an analysis is given of the frequency response characteristics of the temperature-control system 
used in the subcritical suspension reactor at Arnhem (Netherlands), previously described in 


this journal. 


Introduction 


Equipping a plant with a system for automatical- 
ly controlling a variable such as temperature, liquid 
level or gas flow amounts essentially to giving it 
negative feedback and amplification. The output 
signal O, of the instrument that measures the 
variable is compared with a reference (input) signal 
Q;, and their difference ¢ is amplified and used to 
control the regulating device ( fig. 1). It is therefore 
not surprising that control theory runs parallel 
with the theory of negative-feedback amplifiers and 
network analysis, developed in the ’twenties for the 
purposes of electronic engineering !). 


' 5425 


Fig. 1. Block diagram of an automatic control system. I plant. 
II instrument which measures the quantity to be controlled. 
The output signal ©, from II is compared in III with the 
reference signal 0;. The difference « between 0, and O; is 
the input signal of the controlling unit IV, the output signal 
of which controls the regulating unit V. The latter in turn acts 
directly on the quantity to be kept constant. 


An analysis of the characteristics of an automatic 
control system is concerned with the closed loop 
(cf. fig. 1) formed by the plant in which the process 
to be controlled takes place, the measuring unit, the 
controlling unit, and the regulating unit. If all the 
elements composing the control loop are linear, 
the problem can be treated exactly with the aid of 
complex variables. Other methods will not be con- 


sidered here. 


*) Research Laboratories, Eindhoven. 
1) This theory will be found in e.g. H. W. Bode, Network 
analysis and feedback amplifier design, Van Nostrand, New 


York 1945. 


When applying this method to the analysis of a 
negative-feedback amplifier, the starting point is to 
take the characteristics of the amplifier without feed- 
back. The same procedure is applied to a control 
system. The loop is regarded as opened at some point 
— e.g. at the dotted line in fig. 1 — and an 
open-loop transfer function is assigned to it. This 
function, as in AC theory, is the ratio of the 
complex quantities which can be formed from a 
sinusoidal input and output signal — both having 
an angular frequency w — by replacing the trigo- 
nometric function by an exponential function with 
an imaginary exponent. The symbols 0; and 0, 
used above for the input and output signals in fact 
denote these complex quantities. 

The transfer function is normally split into a 
constant real factor K and a complex variable 
G(jw). The latter describes the variation with @ 
both of the gain (amplitude response) and of 
the phase shift y between input and output signals. 
Since all elements of the loop are assumed to be 
linear, KG is equal to the product of the trans- 
fer functions K,G,, K,G,, etc. of the individual ele- 
ments. It is therefore immaterial where the loop 
is imagined to be opened. 

In this article we shall consider three examples of 
the way in which a transfer function is found in 
practice from an analysis of the physical data of a 
given process. The first example relates to control- 
ling the level in a tank, a problem frequently encoun- 
tered in the chemical industry; the second concerns 
the control of temperature in a glass-furnace feeder 
which was recently described in this journal ”), and 
the third the system of controlling the suspension 
temperature in the subcritical nuclear reactor at the 


2) See P.M. Cupido, Philips tech. Rev. 22, 311-319, 1960/61 


(No. 9/10). 
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N.V. KEMA laboratories in Arnhem, also recently 
described in this journal *). It will be seen that the 
physical analysis referred to must often rely on ap- 
proximations. 

For the convenience of readers unfamiliar with 
the subject, we shall first examine the KG function 
at somewhat greater length in order to show how 


fundamental it is to control theory. 


The transfer function 


The importance of the KG function in control 
theory becomes immediately evident when we try 
to answer the following two questions: how far 
does the control system reduce a disturbance Os, 
and is the closed loop sufficiently stable? We will 
first consider the case where the disturbance can 
only occur behind the last block in fig. 1. The 
control loop can then be reduced to the diagram 


in fig. 2. In that case 
Os= KGece O3,, J. 42) 


or, since «= 0; — 0, 


(1 + KG)0, = KGO;, + Os, 
hence 


KG 1 


4 v7 eK C 


LY» tetra 
imebawe oy (2) 


10) 


If we only want to know the extent to which a 
sinusoidal disturbance is suppressed, we may put 
0; = 0, giving 
Oo 1 
-=- on Wet d roe Sa) 
ee Ah Se 


The absolute value of this ratio is called the devi- 
ation ratio *). 
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Fig. 2. Block diagram of control loop in which the disturbance 
@; can only occur behind the last unit (ef. fig. 1). 


It will thus be clear that the effect of the disturb- 
ance is not eliminated, but reduced by the factor 
given in eq. (3). The remaining effect is called the 
deviation. (In theory a zero deviation can be achieved 
by inserting an integrator in the loop, for which 
element G approaches infinity as w approaches zero.) 


*) See B. L. A. van der Schee and M. van Tol, Philips tech. 
Rev. 21, 121-133, 1959/60. 


*) J. M. L. Janssen, Trans. Amer. Soc. Mech. Eng 
mete c. Mech. Engrs. (ASME) 
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As regards the stability of the system the magni- 
tude of @, is unimportant, and (2) reduces to 
tee () 

QO; 1+ KC 

The significance of this equation will be made clear 


presently. 

We now turn to the case where the disturbance 
occurs not at the end of the loop but at a point 
somewhere in the middle (fig. 3). We represent the 
transfer function of the blocks in front of that point 


5427. 


Fig. 3. Control loop where the disturbance signal 0, occurs 
somewhere in the middle. The transfer function of the units 
in front of this point is K,G,, that of the others K,G,. 


by K,G,, that of the others by K,G,. Writing 
K,G,K,G, as KG, the output signal in this case is 


Oo =. Kea KxG, On, ae) 


which, in the same way as above, leads to 


KG KiGs 


Ch —— (OF 
egy Reda iE 


O32 (6) 
From this relationship we again first derive the 
deviation ratio, i.e. the quotient of the values of 


the output signal with and without control action 


for the case where 0; = 0. We find 
Oo il 
: =— —- ‘ree MC) 
0;K.G, 1+ KG 


For the stability of the system we again assume 
QO; = 0 in equation (6), which once more yields 


Oy The 
° = LA ee) 
6; 1 KG 


Apart from the fact, which is trivial, that equations 
(4) and (8) are identical, we see that the right-hand 
side of eq. (7) is identical with that of (3). The 
general rule is that to analyse both the stability of the 
system and the factor by which disturbances are 
reduced it is necessary to have the transfer function of 
the whole assemblage of units in series, i.e. the KG 
function of the open loop. 

The behaviour of an open control loop can also 
be represented graphically, either by a polar plot 
of the KG function in the complex plane (Nyquist 
diagram) or by first deriving from KG the separate 
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amplitude and phase characteristics and plotting 
these in some convenient way. The Nyquist diagram 
is useful for finding quickly the absolute value P 
of the function 1 + KG at a particular frequency. 
This is done by joining the relevant point on the 
curve to the point (—1,0) on the axis, as illustrated 
in fig. 4. 

The frequency response curves are often plotted 
on a logarithmic scale: the amplitude characteristic 
in logarithmic coordinates and the phase charac- 
teristic in semilogarithmic coordinates (Bode dia- 
gram). This method offers advantages when the 
maximum permissible gain is to be determined. 

As far as stability is concerned it will be sufficient 
here to note that the frequency response curves 
should be such that the gain |KG(jw)| is less than 
unity at the frequency w,. where g is —180°, ice. 
where the negative feedback in the loop has changed 
to positive feedback (regeneration). In the Nyquist 
diagram this means that the curve must not en- 


close the point (—1,0) 5). 


Fig. 4. Nyquist diagram of the function KG = K(1 + jwrt) +. 
The distance from a point on the curve to the origin gives the 
absolute value of KG at the relevant frequency, i.e. the gain. 
Disturbances are reduced by the factor P, found from the 
distance of a point on the curve to the point (—1,0). The phase 
shift @ between output and input signal runs here from 0 
to —90°. 


In addition to Bode and Nyquist diagrams the graphical 
methods shown in fig. 5 are sometimes used. The plot in fig. 5a, 
called a Nichols chart, is in effect a Nyquist diagram drawn 
in semilogarithmic rectangular coordinates. In fig. 5b, called 
the inverse Nyquist diagram, the reciprocal of KG is drawn 
in the complex plane. This method is useful for quickly deter- 
mining the response of the closed loop, since it follows from 


(4) that 0,/O. = (1 + KG)/KG = KG + 1. 


5) For complicated control systems the stability criterion 
as formulated here is oversimplified, but it is adequate for 
the relatively simple systems discussed in this article. The 
stability problem is dealt with at some length e.g. by 
G. S. Brown and D. P. Campbell, Principles of servome- 
chanisms, Wiley, New York 1948, and by G. J. Thaler 
and R. G. Brown, Servomechanism analysis, McGraw-Hill, 
New York 1953. 

The relation between stability and gain will be dealt with 
in the next number of this journal. 
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Fig. 5. a) Relation between log |KG| and ¢ for the transfer 
function KG = K(1 + jor) (cf. fig. 4) plotted in rectangular 
coordinates (Nichols chart). 

b) Graph of the function 1/KG in the complex plane, likewise 
for the case KG = K(1 + jor). 


Example of the determination of the transfer func- 
tion of a single element 


As an example of the method used to determine 
the transfer function of an element we take an open 
tank into which water is run from a tap (rate of 
flow I,) and at the same time run off through an 
outlet (rate of flow I,). The tank is assumed to be 
of uniform width, area of base F’, so that the volume 
of water V in the tank is proportional to the height 
H to which the tank is filled (fig. 6). It is further 
assumed that I, is proportional to H, which we 
may write I, = cH. If I, is constant (value I)), 


V5 
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Fig. 6. For calculating the transfer function of a tank filled 
with water from a tap (flow rate J,, independent variable) 
and emptied through an outlet where the flow rate I, is pro- 
portional to the height H. 
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the value H, which H has in the steady state is 
such that J, is equal to I), and therefore Hy = 
I,/c. We now want to ascertain how H will fluc- 
tuate around the value Hy) when J, is not constant 
but shows a sinusoidal ripple of amplitude 1, and 


angular frequency ©: 
I, = Io ty sin wt. 
The problem is defined by the differential equation 
dV dH 
Se eae ee 
dt dt 


(ee 


or, substituting H, + h for H and given I, = cH), 


dh A 
F— + ch=i,sinot. fp ee(9) 
dt 


We can readily solve an equation of this kind by 
finding a solution of the form h=he™ for the 
equation 


(10) 
where i, = i, e!”’. Substituting the expression for 
h in eq. (10), we obtain 

Fjoh +ch= ae 
or 


es a bale 


Fj SEG ‘ij es jor’ 


(11) 


where t represents the quotient F'/c, which is easily 
seen to have the dimension of time. Since finding 
the KG function consists in finding the ratio of the 
complex quantities corresponding to the input and 
output signals, it follows from equation (11) that 


1/e 
KG = ———_., 
1 + jor 
whence 
1 
K= Ve and €=——— =. (12) 
1 = jor 


(The factor 1/c in the numerator of (11) is not a 
gain factor but a dimensional factor, due to our 
choice of fh as the starting quantity instead of the 
flow I,. If we had started with I,, K would have 
been unity; there is no amplification of the input 
signal.) 

The form of the G function found here is frequent- 
ly encountered in practice. The operation of most 
circuit elements used in control systems can be 
defined by an arrangement of blocks in series — as 
many blocks as there are time constants — having 
a transfer function of the form KL jon). 
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Table I. Electrical circuits whose transfer function can be 
expressed in the form K(1 4 jr)’. In the third example the 
input and output quantities have different dimensions, result- 
ing in a dimensional factor R. The gain factor is in each case 
equal to unity. 


Input 


and 
Circuit output Transfer function %E 
quan- 
tities 
R 1 
—______ = ________._ | IR 
Vos R+job 1-+joL/R / 
R L/joc il 
A Vy, ty We = RC 
’ ia 0 Wateae I/jjoC-+R 1+ joRC 
Lf Vy 
c Rela 
WR oe 1 + joRC 


Some familiar electrical circuits that have this trans- 
fer function are shown in Table I. The last one may 
be regarded as the equivalent electrical circuit of 
the tank in our example. 


Control of the glass temperature in a tank-furnace 
feeder 


As our second example we shall consider a system 
for controlling the temperature of molten glass 
in the feeder from a tank furnace to a glass-working 
machine. In the article cited *) it was explained that 
it is important to ensure that the conditions under 
which glass is produced are kept as uniform as 
possible. The first prerequisite is to keep the tem- 
perature of the glass in the melting tank carefully 
constant. Since this cannot be done directly, the 
various quantities which affect the glass tempera- 
ture are stabilized separately. These quantities are 
the temperature of the fuel oil for the burners, the 
pressure of the atomizing air, the flow rates of oil 
and air to the burners (or possibly the flow rate of 
the air and the ratio of oil flow to air flow), the pres- 
sure of the combustion gases above the glass, and 
the level in the tank. From the melting tank the 
glass flows through the “working end” of the fur- 
nace into one or more feeders for the glass-working 
machines. Since the flow of glass to the machines 
must be constant, the temperature in the feeder is 
separately controlled; this control must be very 
rigorous, because the viscosity of the glass, which 
has a great effect on its rate of flow, is highly 
temperature-dependent. We shall now examine the 
control system used for this purpose. 

The temperature of the glass in the feeder is con- 
trolled with the aid of an electrical thermometer, 
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which is immersed in the molten glass and whose 
output signal is used to regulate the burners — gas 
burners in this case ( fig. 7); we shall not discuss the 
practical details of this control here. A block dia- 
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Big. 7. Heating of molten glass in a glass-furnace feeder. 
I glass, 2 wall of feeder. 3 roof. 4 gas burner 


gram of the control system is shown in fig. 8. Block I 
is the feeder with molten glass, which receives a spe- 
cific flow of heat and acquires the temperature Oj. 
Block II is the thermometer, whose output signal 
Oy is compared with the desired value 0;. The 
difference ¢ is amplified in block ITT — we assume 
provisionally that the amplification A is linear and 
independent of frequency. The output from III 
controls the valve IV which regulates the gas supply 
and thus indirectly the flow of heat to the gas. Let 
us now try to find the transfer function of this 
system. 

We shall argue presently that the feeder with 
molten glass, like the water tank discussed, behaves 
in principle like the electrical circuit of a resistance 
and capacitance in parallel (see example 3 in Table I). 
The current I represents the flow of heat from the 
burners to the glass, the capacitance C the heat 
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Because of the very high values assumed by the 
heat capacity of the glass and the thermal resist- 
ance, T, is of the order of magnitude of half an hour. 

In order to justify the assumption that this case may 
in fact be identified, as a first approximation, with 
the third example in Table I, we shall look at what 
happens in the feeder (see fig. 7). The heat deliv- 
ered by the burners is partly transferred to the 
glass and partly to the walls and roof of the feeder. 
The ratio between the two heat flows will be roughly 
independent of the magnitude of the total heat 
flow. The way in which the heat is transported from 
the flame to the glass need not be considered here. 
All we are interested in is the fact that there is a 
flow of heat to the glass that depends on the flow 
of gas to the burners, and which responds almost 
without any lag to a change in the gas flow. 

It remains to be shown that the process of heat 
removal may be regarded as the flow of current 
through a resistance as a result of a potential dif- 
ference Vo which represents the temperature dif- 
ference AT between the glass and the surroundings 
of the furnace. The heat is transported partly by 
radiation and partly by conduction and convection: 
the glass radiates heat to the roof and walls, heat 
is conducted through the glass and walls to the 
atmosphere, and finally heat is lost from the out- 
side of the feeder by radiation, conduction and 
convection. The dependence of the heat flow on 
the temperature difference AT can, of course, be 
expressed in a power series. This series, owing in 
particular to the presence of a radiation component, 
will also contain higher powers of AT. It follows, 
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Fig. 8. Block diagram (simplified) of system for controlling the temperature 
of the molten glass flowing through a feeder from the glass furnace to the 
glass-working machine. I feeder. JJ thermometer. IJJ linear amplifier. 


IV regulating valve. 


capacity of the glass, R the thermal resistance over 
which the heat is dissipated, and V_ the temperature 
to be kept constant. Disturbances that arise may 
be due to fluctuations in I, C or R. Fluctuations in 
C may be due to variations in the level of glass, fluc- 
tuations in R to variations in the chimney draught. 
We shall confine ourselves here to variations in I. 
Putting RC = 1, the frequency-dependent part 
of the transfer function for block I is given by 


1 


ee, 13 
1+ jor, o?) 


Gr 


then, that the thermal resistance, i.e. the derivative 
of the heat flow with respect to AT, is not inde- 
pendent of AT and in an equivalent electrical circuit 
may only be represented by a pure resistance 
provided AT is practically constant. In our case 
this condition is satisfied, the temperature varia- 
tions of the molten glass in the feeder being of the 
order of only 1 °C. 

Not so much need be said about the other units 
of the loop. The thermometer, block IJ, can be 
treated as a capacitor (heat capacity) ¢ which is 
charged by a current supplied by a voltage source 
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via a resistance r. Here the input quantity is the 
voltage which represents the temperature of the 
glass tank (Table I, example 2). The transfer func- 
tion for G,,, like that for G;, therefore again has the 
form (1 + jwt,)~!, but the value of the time con- 
stant tT, (= re) is very much smaller (e.g. half a 
minute). 

The transfer function of block IV, the regulating 
valve, must define the relationship between the 
output signal from III and the flow of heat to the 
tank. There are two steps involved. The signal from 
III controls the valve in the gas pipe. The relation- 
ship between this signal and the corresponding 
change in the flow of gas will not, in general, be 
exactly linear. Secondly, the flow of gas to the 
burner determines the flow of heat to the glass; 
the relationship between these two quantities, 
which is partly governed by the size of the luminous 
cone in the flame °), will not always be linear either. 
For small variations, however, the behaviour of a 
non-linear element like block [V may permissibly 
be regarded as linear. Since a valve usually shows 
a certain lag, which may be characterized by one 
time constant, we have G,,=(1 + jt). Here 
again, T;y is small compared with 17,. 

Let us now briefly review the simplifications we 
have adopted. Firstly, of course, there were the 
approximations introduced to simplify the for- 
mulae for G; and Gy. Further, the block diagram 
as such (fig. 8) contains various approximations. 
In the first place we have disregarded the fact that, 
owing to the heat capacity of the burner, the heat 
flow to the glass does not respond instantaneously 
to a change in the gas supply. To take this into 
account we must add the factor (1 + jmty)~! to 
the transfer function of the whole loop. It is also a 
simplification to represent the glass by a single 
capacitance and to neglect the heat capacity of the 
wall. Finally, the temperature distribution of the 
glass is certainly not entirely uniform. There is no 
space to allow for all these factors here. 

As a first approximation, however, our simplified 
treatment of the system’s behaviour is justified. 
We may even go a step further and, since 1,, and 
Try XK Tj, equate G,, and Gy, to unity. The block 
diagram then reduces to that of fig. 9, and the trans- 
fer function of the open loop is 

KG aon 
1+ jor, 


where K — A, and G = (l + jot;)>. 


(14) 


°) See page 313 of reference 2). 
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As far as the stability of the system is concerned, 
however, the presence of the small time constants 
can not be neglected, for the variations of the con- 
trolled quantity are required to be very small, and 
therefore K must be given a very high value (see 
eq. (3)). The extent to which that is possible without 
endangering the stability of the system is in fact 
governed by the value of these small time con- 
stants ”). 
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Fig. 9. Simplified form of block diagram of fig. 8, obtained by 
equating with unity the G function of all blocks for which 
t <tr. This approximation is permissible only if the gain 
factor A of III is not too high. 


Control of the suspension temperature in the KEMA 
subcritical suspension reactor at Arnhem 


Before finding the open-loop transfer function for 
the subcritical suspension reactor at Arnhem, we 
shall briefly recapitulate the operation of the control 
system, which has already been described in this 
journal *), A suspension of uranium oxide is pumped 
at a constant rate around the reactor circuit (see 
fig. 10). The pump P supplies 3-7 kW of energy 
to the fluid, and therefore some cooling is necessary 
even at the highest operating temperature. Sta- 
bilization of the temperature against variations 
in the power supplied by the pump and withdrawn 
by the cooler is effected by applying somewhat ex- 
cessive cooling and by compensating the excess 
with the aid of a variable heater element V, which 


Fig. 10. Schematic representation of reactor circuit (left) and 
cooling circuits of the KEMA subcritical suspension reactor 
at Arnhem. R reactor vessel. P pump which circulates the 
suspension of uranium oxide around the reactor circuit in the 
direction of the arrows. T resistance thermometer. C control 
unit. V heater (regulating circuit). W, first heat exchanger. 
I primary cooling circuit. W, second heat exchanger. 2 sec- 
ondary cooling circuit. K cooling-water valve. 


") To be discussed in the forthcoming article mentioned in 
footnote °), é 
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surrenders its heat directly to the suspension. The 
heater element is controlled by a thermometer T, 
which measures the temperature of the suspension. 
At equilibrium, the heater gives up about half its 
maximum power (5 kW), making it possible to 
offset variations of more than 2 kW in the power 
to be dissipated. The heater is simply a length of 
metal piping through which an electric current is 
passed. The temperature is measured with a 
Wheatstone-bridge resistance thermometer which 
delivers a voltage proportional to the difference 
between the measured and the desired temperature. 

A block diagram of the control system is shown 
in fig. 11. The first block (A) again represents a 
linear amplifier. The transfer functions of blocks I 
and ii together define the behaviour of the heater 
(i.e. the way in which the temperature 0; of the 
effluent suspension reacts to a variation in the power 
supplied to the heater), that of the last block (ITI) 
the behaviour of the thermometer itself. The dis- 
turbance @, is a change in the temperature of the 
suspension, which is not detected until liquid of 
the wrong temperature reaches the thermometer. 
A fluctuation in the mains voltage must also be 
introduced at this point. 
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Fig. 11. Block diagram of temperature-control system in the 
Arnhem subcritical suspension reactor. Blocks I and IJ to- 
gether represent the heater, block III the thermometer. 


It may be asked whether fig. 11 does in fact represent the 
control loop, for there is another closed loop involved, viz. 
that of the circulating suspension. In order to show that the 
latter loop may be disregarded for our purposes, we recall 
that the control system has been designed to offset variations 
in the power delivered to the suspension by the pump and 
removed by the cooler *). The control equipment need only 
ensure that the fluid leaves the heater at the correct tempera- 
ture. It is immaterial whether or not the fluid itself flows in an 
open or closed circuit. 

Nor are our considerations affected by disregarding the 
influence of heat losses in other parts of the reactor circuit. 
Any such variation can be treated as a disturbance signal en- 
tering the loop at the same point as 0 , and can thus be 
regarded as a component of the latter. 


To determine the transfer function of the heater 
it is convenient to divide the heating process into 
two stages, represented in fig. 11 by blocks Tand II. 
The transfer function K,G, defines the response of 
the temperature of the wall of the heater to a va- 
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riation in the supplied power; the function K,,G,, 
defines the response of the temperature of the ef- 
fluent to a change in the temperature of the wall 
of the heater. In this connection we make two 
closely related simplifying assumptions. Firstly, 
we assume that the temperature of the heater pipe 
is uniform at all points, so that we may in fact refer 
to the temperature of the pipe. Secondly we assume 
that variations in the temperature of the effluent 
are negligible compared with the temperature va- 
riations shown by the pipe wall; in other words, we 
regard the heat capacity of the fluid as extremely 
high. In fact the suspension circulates fast enough 
for these assumptions to be correct. Otherwise, our 
division of the heater into two blocks would not be 
permissible, unless feedback were applied from 
block IJ to block I. 

The transfer function K,G; can most quickly be 
found from the equivalent electrical circuit shown 
in fig. 12a. The current I represents the power 
supplied, C, the heat capacity of the pipe, R, the 
thermal resistance at the boundary between pipe 
and suspension, C, the heat capacity of the suspen- 
sion, V, the temperature of the pipe and V, the 
temperature of the fluid. Resistances R, and Ry 
represent heat losses due to conduction and con- 
vection in the pipe and fluid mass in the rest of the 
reactor circuit. Measurements show that C, ~ 
500 J/°C, C, ~ 400 C,, R, ~ 25x 10-3 °C/watt and 
R, ~ 2.2 10-3 °C/watt. We see that C, is in fact 
much larger than C,. We can thus regard C, as 
short-circuiting R,, in which case the diagram in 


Vp 
C2 = Rs 
pi Yj ls Yj 
Ry ay) Ro Cy 
iz b bd Cc 5431 


Fig. 12. a) Equivalent electrical circuit of heater operation. 
As C, > C,, this can be simplified to 6) or, by substituting 
R, for R, and R, in parallel, to c). 
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fig. 12a reduces to that in fig. 12b, and the latter 
may in turn be simplified to fig. 12¢ by replacing 
R, and R, by a single resistance Ro. 
The transfer function of the first block is there- 
fore equal to (see Table I, example 3): 
Va Hie pees 
I 1+ jor, 
ie, Kp = Ri = 2 °C/KW and Gio) = > jou)=- 
The time constant 7, is equal to RoC, (one second). 
The transfer function K,,G,, (the change in the 
temperature of the fluid due to the change in the 
heater temperature) is found by determining the 


(15) 


way in which the temperature of the effluent sus- 
pension fluctuates when the temperature of the 
pipe varies sinusoidally around the equilibrium 
value T, with an amplitude T, which is small 
compared to the average temperature difference 
between the pipe wall and the suspension. (Since 
the flow is highly turbulent, the temperature of the 
fluid is virtually identical in all points of a cross- 
section.) 

If we consider a “disc” of suspension of thickness 
dx, moving at a constant velocity v through the 
pipe, we see that the disc must constantly encounter 
a different wall temperature in its passage through 
the pipe. The difference between this temperature 
and the equilibrium temperature of the wall is 
T, sin w(t) + x/v). Here t) is the time at which the 
disc entered the pipe, and x the distance which the 
disc has travelled in the pipe at the given moment. 
We shall now consider how much extra heat is 
supplied to the disc in its passage through the pipe 
as a result of this temperature difference. In the 
time dx/v during which the disc remains between 
x and x + dx it takes up an additional quantity of 
heat which gives rise to a temperature change 
proportional to 


(16) 


; 55 dx 
T, sin (ty ae =) SCE 
v v 
The integral of this expression over the whole length 
! of the heater is proportional to the difference 
between the actual temperature which the disc 
has acquired by the time it reaches the end of the 
pipe and the temperature it would have acquired 
if the wall temperature had permanently been equal 
to T). This difference is thus proportional to 
l 


T. si x\ dx 
1 Sin w(t) + -) - (17) 
; v/) v 
0 
which reduces to: 
: SNOT, —. 
MGyy Ll STOTT (18) 
OO Tat 
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where 1,, = I/2v (about 0.15 second). If we now 
consider the entire flow, i.e. a continuous series of 
discs for which the value of t,) is successively larger 
by an amount dx/v, we may conclude that the varia- 
tion of the temperature difference with time is given, 
apart from a constant factor, by expression (18) if 
ty is replaced by (¢ — l/v), which is equal to (t—27%,). 
Expression (18) then becomes: 

sIn@ Ty 


2t, T,— 
OT 


(19) 


sin w(t — Ty) - 


The corresponding complex quantity (we shall 
henceforth disregard the constant factor 27,;) is: 


sinw T : 
EOE rn)s 


ata O40) 


OT 
To find G,(jw) we divide (20) by the complex quan- 
tity corresponding to the input signal, T,e), 
giving: 


GAljo) =e (21) 


OTT 


It can be seen from eq. (21) that the time 1; (which 
is known as the distance velocity lag) is of the nature 
of a transit time. This transit time causes a phase 
shift g = —wt, proportional to o ( fig. 13). Further, 
there is a frequency-dependent amplitude factor 
(sin @T,;)/@t Which is equal to unity at very low 
frequencies and is zero when wt, = 2. In the latter 
case each disc traverses the heater in exactly one 
period, and thus takes up just as much energy in 
the one half as it gives up in the other half. 

Block III in the diagram defines the behaviour 
of the resistance wire of the thermometer. This be- 
haviour can be characterized to a good approxima- 
tion by a single time constant 71; (about 0.10 sec), 
so that 


1 
Gyy(jo) = —— : Ree Kia 
mJ ) ag ee (22) 
——— (a) 
5 10 15 20 25 rad/sec 


Fig. 13. Phase characteristic of the dissipation of heat by the 
heater in the circulating suspension (cf. eq. (21)). Because of the 
distance velocity lag ty, there is a phase shift proportional to 


. At the frequency » = m/t1x (© 25 rad/sec) the phase shift 
p is —2, 
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The transfer function of the open control loop 
as a whole, which is equal to the product 4 K,G,(ja) x 
KyGy(jo) KyyGyy(jo), is thus: 

C) 1 
— = AK" Sd Site = 
E 1-- jw, 


SLO) tre i} 
—_—_ e _JO@TII SZ 
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(Since the thermometer is not exactly at the end 
of the heater but somewhat further up in the pipe 
line, we should introduce into eq. (23) a second 
distance velocity lag, equal to the time t,y which 
the fluid needs to cover the relevant distance. How- 
ever, since Ty <tTy, this factor may be neglected.) 

The frequency response curves derived from this 
transfer function are shown in fig. 14, and the Ny- 
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Fig. 14. Frequency response curves of the open control loop. 
Top: amplitude characteristic. Bottom: phase characteristic. 
At very low frequencies the gain is roughly 7. At the value 
of w where the phase shift ¢p is 180°, the gain is less than unity. 


quist diagram in fig. 15. In fig. 16 the factor P, 
the reciprocal of the deviation ratio, is plotted as a 
function of o. 

A word should be added about the fact that the 
system is unable to correct high-frequency disturb- 
ances (see fig. 16). A sudden change in the valve 


APPLICATION OF LINEAR CONTROL THEORY 117 


1 


Fig. 15. Nyquist diagram pertaining to fig. 14. 


controlling the cooling water in the secondary 
cooling circuit will not result in an immediate 
change in the temperature of the liquid flowing 
over the thermometer, nor will a sudden change in 
the mains voltage, i.e. in the power fed back to 
the heater (cf. fig. 12). Rapid variations of this kind 
thus give rise to relatively slow disturbances in the 
temperature of the suspension — which will, of 
course, be corrected — and are therefore unable to 
do any harm. 

In conclusion, some general remarks. The above 
examples have shown that it is very often not 
possible to calculate a loop transfer function exactly, 
and that approximations have to be made. In any 
case, exactly linear systems are never in fact en- 
countered, so that our assumption of linearity is in 
itself an approximation. In cases where it is not 
possible to calculate the KG function at all, the 
frequency response curves can in principle be found 
from measurements. A sinusoidal disturbance is 
deliberately introduced, and the resultant variation 
of the output signal (OQ, in fig. 1) is then recorded. 
This method cannot of course be used where the 
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Fig. 16. The factor P (cf. fig. 4) as a function of w, derived 
from fig. 15. 
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system is non-linear, or where the control loop can- 
not be opened. Considerable practical difficulties 
arise where long time constants are involved, be- 
cause one then has to wait a very long time in each 
measurement before the switching transient has 
ended. In the latter case the object can sometimes 
be achieved by introducing a step-function disturb- 
ance and seeing how the system then behaves. 
From the step-function response thus obtained one 
can derive the frequency response curves, and from 
these find the values of the various time constants. 
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Summary. The theory of an automatic contro Isystem is formally 
similar to that of a negative-feedback amplifier. The behaviour 
of a feedback system, provided all its elements are linear, can 
be derived from the open-loop transfer function, KG(jq). 
This function is the quotient of the complex quantities corre- 
sponding to a sinusoidal input signal (angular frequency ©) 
and the resultant sinusoidal output signal. K is a constant 
factor. Disturbances are reduced by the factor |1 + KG]. 
The system is stable when the curve representing KG in the 
complex plane does not enclose the point (—1,0). The KG 
functions are calculated for 1) a tank with water flowing 
in and out, 2) the temperature control system for a_glass- 
furnace feeder, and 3) the temperature control system for the 
subcritical suspension reactor at Arnhem. The elements of a 
control loop frequently have a transfer function of the 
form KG(jw) = K(1 4+ jot). 


A SIMPLE ANALOGUE COMPUTER FOR DETERMINING 
THE COLOUR POINT OF A LIGHT SOURCE 


by B. van der WAAL *). 


535.651.1 


Determining the “colour point’ of a light source calls for a series of simple but laborious 


computations. The computer described below performs this work rapidly and automatically ; 


the CIE-distribution coefficients, which define the spectral sensitivity of the eye of a 


“standard observer’’, are built into the machine. 


One of the principal characteristics of a light source 
is the colour of the light it emits. The colour can be 
objectively specified!) with the aid of three quanti- 
ties X, Y and Z, defined as 


X = J S(A) x(A) da, 


VS eS) Cl) die ee eet} 
i, MOV AM OV 


In these equations S(A) is the relative spectral energy 
distribution of the radiated light; x(/), y(A) and 
2(A) are the CIE-distribution coefficients, which de- 
fine the spectral sensitivity of the eye as laid down 
in 1931 for the “colorimetric standard observer” by 
the International Commission on Illumination (CIE). 
The coefficients are so chosen that the values of y(A) 
are identical with those of the relative luminous 
efficiency V(A). The maximum values of these coef- 


*) Lighting Division, Eindhoven. 

1) See W. de Groot and A. A. Kruithof, The colour triangle, 
Philips tech. Rev. 12, 137-144, 1950/51; F. W. de Vrijer, 
Fundamentals of colour television, Philips tech. Rev. 19, 
86-97, 1957/58; P. J. Bouma, Physical aspects of colour, 
Philips Technical Library, 1947. 

The notation used in these articles differs somewhat from 
current usage, which conforms more or less with the 
official recommendations of the International Commission 
on Illumination; see “International Lighting Vocabulary”’ 
publication CIE 1-1-1957. i ; 


ficients occur at wavelengths in the red, green and 
blue parts of the spectrum, respectively (see fig. 1). 

It is however usually more convenient to specify 
the colour of a light source by the ratios of X, Y and 
Z to their sum: 


As the sum x + y + z is 1, two of these values are 
of course sufficient to specify a colour. Any colour 
can thus be represented by a point in an x-y dia- 
gram, as in the CIE chromaticity diagram or colour 
triangle (fig. 2). The colour points of monochromatic 
light lie on a horseshoe-shaped curve called the spec- 
trum locus. The straight line joining the ends of the 
spectrum locus, the “purple line’’, represents the pure 
(saturated) purples, mixtures of red and violet radia- 
tions. In the centre lies the “white point” FE, where x, 
y and z are equal. The other points may be regarded as 
colour points of a mixture of white light with a spec- 
tral colour or with a colour on the purple line. 


2 Oe ne tonal 
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700 nm 
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Fig. 1. The CIE-distribution coefficients x(A), y(A) and z(A), 
giving the equal-energy distribution curves in the CIE 
system. 


(A chromaticity chart showing the actual colours 
corresponding to various points can be found 
in the first article referred to in note }). 


The determination of the colour point of a light 
source 


The method of finding the colour point of a given 
light source follows directly from equations (1) and 
(2). The value of S(A) for the light is measured with 


1,00 
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Fig. 2. The CIE chromaticity diagram in x-y coordinates. On 
the spectrum locus (curve of the spectral colours) pane 
sponding wavelengths are given in nm. (lial Omen) 
E is the “white point’’. The curve inside the spectrum locus 
joins the colour points of a black-body radiator at different 


temperatures. 
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a spectrophotometer at a series of wavelengths rising 
in steps of 10 nm from 385 to 735 nm. (lnm = 10-9 m.) 
To find X the values thus measured are multiplied by 
the values of x(A) corresponding to the same wave- 
lengths and the products are then added. (One should 
really multiply by the width of the wavelength inter- 
val used in the measurement, but as we are only con- 
cerned with the ratios of X, Y and Z, we have omit- 
ted to do so here.) The values of Y and Z are cal- 
culated in the same way, after which x and y are 
determined by two division operations from equa- 
tion (2). The computation of a single colour point 
thus involves about 100 multiplications, followed 
by three additions and two divisions. 

A simple analogue computer has been designed 
for performing these operations automatically. The 
principle is illustrated in fig. 3. For every wavelength 
An at which S(A) is measured there is a circuit con- 
sisting of three resistors Ryp, Ryn and Rn, whose 
respective values are inversely proportional to 


x(A), y(A) and 2A) at that wavelength. A voltage 


Fig. 3. Principle of the computer circuit for determining the 
colour point of a light source. 


derived from an accumulator B, via a potentiometer 


Rp 


age being proportional to the measured value of 


nis applied to each resistor, the value of the volt- 


S(A) (for this purpose the meter MV, used as a volt- 
meter, can be connected to each potentiometer in 
turn). The currents through the resistors are then 
proportional to S(An)*(An), S(An)y(An) and S(An)z(An)> 
respectively. The sum of the currents through all 


resistors Ry» is 
=S(An) (An) 


and is therefore proportional to X. The sums of the 
currents through Ry, and R;, are similarly propor- 
tional to Y and Z respectively, the constant of 
proportionality being the same in all cases. 
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In measurements on tubular fluorescent lamps 
the energy at the wavelengths of the five mercury 
lines is determined in addition to S(A) at the wave- 
lengths mentioned. The computer contains similar 
circuits for these wavelengths too, giving a total 
of 41. 

When all 41 potentiometers have been set to the 
appropriate values, the three currents are deter- 
mined by a null method, using a circuit consisting 
of a galvanometer G, a battery B,, the meter M 
(now used as an ammeter) and a variable resistance 
R,. The switch S, is closed, S, opened, and S, 
turned first to position X. The current through G 
is now made zero by adjusting R,; X can then be 
read off from meter M. The values of Y and Z are 
found in the same way, with S, in positions Y and Z 
respectively. 

Switch S, was given the form shown in fig. 3 to 
ensure that all resistors pass current both during the 
adjustment of the voltages and the measurement 
of the currents. If current were to flow through only 
one of the three sets of resistors (in which case only 
the top half of S,; would be needed) the currents 
through the potentiometers would change when S, 
was switched over, and the voltages would have to 
be reset. With this arrangement, however, one set- 
ting serves for all three measurements. A photo- 
graph of the computer can be seen in fig. 4. 


Other uses for the computer 


The computer can be used for many other cal- 
culations that involve multiplication operations 
with the CIE-distribution coefficients or simply 
For 
example, to find the distribution of luminous flux 


over various wavelength regions, the measured 


with the relative luminous efficiencies V(A). 


values of the relative spectral energy distribution 
S(A) are multiplied by the values of V(A) for the 
wavelengths concerned, and the products for the 
various wavelengths are added. 

The computer is also useful for “normalizing” the 
relative spectral energy distribution S(A). This is 
done when it is necessary for irradiation purposes 
(e.g. the irradiation of plants) to know the absolute 
value of the radiant power delivered to a surface by 
a light source in a particular range of wavelengths. 
A spectral energy distribution is then needed which 
gives the absolute instead of the relative values for 
the surface in question. 

These values might be determined by the follow- 
ing procedure. The total radiant power P incident 
on the surface is measured and divided by [S(A)d4, 
an integral proportional to the total radiant power 


received by the photometer in the measurement of 
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Fig. 4. Photograph of the colour-point computer. The sloping 
panel contains the potentiometer controls for presetting the 
measured energy values. Above left, the meter M shown 
in fig. 3, and right, the galvanometer G. The bottom row 
contains, among other things, the controls of the variable 
resistors used in the determination of X, Y and Z. 


S(A). Multiplying the values of S(A) by this factor 

P/{S(A)dA gives the absolute distribution required: 
S(A) 

fS(A)da 

S(A)/[S(A)dA is called the normalized distribution. 


To find from this the absolute distribution on an 
arbitrary surface the values are multiplied by the 


C(O) 


measured total radiant power through that surface. 

The actual procedure is somewhat different. For 
practical reasons it is easier to measure the luminous 
flux on a surface than the radiant power. The nor- 
malized energy distribution is then found by divi- 


ding S(A) by 
KSA) VG)\dis =.) een) 


(This integral is proportional to the total luminous 
flux received by the photometer in measuring S(A); 
K is the luminous efficiency of radiation, which 
is about 680 lumen/W.) The normalized energy 
distribution found in this way is expressed in 
uW/lumen per 10 nm. The absolute distribution 
is found from it by multiplying by the measured 
luminous flux on the surface. 

Since V(A) and y(A) are identical, the integral 
in eq. (4) is equal to the integral Y in eq. (1), and can 
thus be determined with the computer in the same 
manner as the colour point. Now, however, the 
absolute value of Y must be known, and therefore 
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the computer is calibrated by determining Y with 
it for an energy distribution in which this quantity 
is known. For this purpose we used the equi-energy 
spectrum, for which Y can be easily computed. 


Summary. After a brief explanation of the CIE chromaticity 
chart, a simple analogue computer is described with which 
the colour point of a light source can be determined from the 
relative spectral energy distribution S(A), measured at a series 
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of wavelengths rising in steps of 10 nm. These values are 
multiplied by the corresponding CIE-distribution coefficients 
x(A), y(A) and 2(A). The products S(A)x(A), S(A)y(A) and S(A)z(A) 
are all added to find the values of X, Y and Z which define 
the chromaticity of the light. For every wavelength at which 
S(A) is measured, the computer contains three resistors whose 
values are proportional to 1/x(A), 1/y(A) and 1/z(A). A voltage 
proportional to the measured value of S(A) is applied to 
these resistors. The currents then flowing through them are 
proportional to the products mentioned, and are added in 
the computer. Other operations which can be carried out 
with the aid of the computer include the normalizing of a 
relative spectral energy distribution. 
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ANALYSIS OF RESIDUAL GASES IN TELEVISION PICTURE TUBES 
WITH THE AID OF THE OMEGATRON 


by J. van der WAAL *) and J. C. FRANCKEN *). 


621.385.832.032.94:621.384.83 


The omegatron, in a new design recently described in this journal, is a particularly useful 


mass spectrometer for the quantitative analysis of residual gases. Analyses of this kind have 


provided valuable information on processes inside a television picture tube during manufacture 


and in operation. 


In a properly functioning television picture tube 
the total residual gas pressure should be less than 
10-> torr (1 torr = 1 mm Hg). The composition 
of the residual gas is also an important considera- 
tion, for one kind of gas can be more harmful than 
another. 

The life of the cathode in a picture tube is par- 
ticularly dependent on the types of gas present. 
The composition of the residual gas in an evacuated 
and sealed-off tube changes continuously, owing to 
the desorption and adsorption of gases by the vari- 
ous components. The composition of this residual 
gas can be determined with a mass spectrometer. 

It is important that the spectrometer used should 
have a small volume and that it should release little 
gas compared with the picture tube. A mass spectro- 
meter found especially suitable for this purpose is 
the omegatron. 

The omegatron had already been used in 1953 
for qualitative analysis of residual gases in sealed- 
off cathode-ray tubes !). As recently mentioned in 
this journal, an omegatron has now been designed 
which is also suitable for quantitative analyses ”). 
Its use for analysing the residual gases in television 
picture tubes will be described in this article. The 
operation of the omegatron was discussed at length 
in the above-mentioned article. For convenience, 
we shall briefly recapitulate this discussion here. 


Principle of the omegatron 


A perspective sketch of the omegatron in its 
simplest form is shown in fig. 1. The narrow beam of 
electrons from the cathode K is parallel to a uni- 
form magnetic field. The collision of electrons with 
gas molecules gives rise to ions. If the latter have a 
velocity component perpendicular to the direction 


*) Electron Tubes Division, Eindhoven. 
1) J. Peper, Philips tech. Rev. 19, 218-220, 1957/58. 


*) A. Klopfer and W. Schmidt, Philips tech. Rev. 22, 195-206, 
1960/61 (No. 6). 


of the magnetic field, they describe circular paths 
perpendicular to the static magnetic field. The radii 
of these paths are given by the equation: 


yeaa ck Nee ere. (Mi) 
ap fe 
where m is the mass of the ion, e its charge, vg the 
component of the velocity perpendicular to the 
direction of the magnetic field, and B the magnetic 
induction. The angular frequency «, of the revolu- 
tion of the ion is 


(he aes Pe es es (P4) 
m 


and the period of revolution is therefore independent 
of the velocity vy. An RF field Ey¢ sin wot is applied 
at right angles to the magnetic field. If the angular 
frequency «, of that field is made equal to the an- 
gular frequency «, of a given kind of ion, the result 


V 


3249 | 


Fig. 1. Perspective sketch of a simplified omegatron. B indi- 
cates the direction of the magnetic induction, Ep, the direction 
of the RF field. The figure illustrates the spiral path described 
by an ion. K cathode, which emits the ionizing beam of elec- 
trons el. A and H are the electrodes to which the RF voltage is 
applied. I ion collector. V connection to amplifier. 
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is a condition of resonance. The resonating ions then 
describe spiral paths of uniformly increasing radius. 
They are caught on a suitably placed collector elec- 
trode, and the ion current, i.e. the number of ions 
collected per unit time, is detected with the aid of a 
sensitive amplifier. 

Ions having a different mass but the same charge 
as the resonating ions, i.e. ions with a different 
value of e/m, are periodically accelerated and dece- 
lerated by the RF field. The radius of the orbits 
described by these ions alternately increases and 
decreases, but the maximum value generally re- 
mains smaller than the distance between the point 
of origin of the ions and the collector. Such ions do 
not therefore contribute to the measured current. 
To discriminate in this way between two kinds of 
ions, their mass difference Am must exceed a spe- 
cific minimum value. The ratio m/Am is called the 
resolution of the omegatron. 

The pressure of each component of the gas mix- 
ture can be found by varying the frequency a, 
of the RF field and measuring the ion current for 
the various mass numbers. If the relationship 
between ion current and pressure is reproducible, 
the partial gas pressure for any mass number can 
be determined from the ion current. This relation- 
ship will be reproducible when all resonating ions 
reach the collector electrode. That is only the case 
when the total gas pressure in the omegatron is less 
than 10~- torr; otherwise ions are lost by collision 
with gas molecules. Further, ions in resonance 
should not strike other electrodes. This can be 
prevented by fitting side plates and appropriately 
choosing the potential distribution inside the ome- 
gatron, as described in the article quoted *). A photo- 
graph of the omegatron used for this investigation 
is shown in fig. 2. 


Experimental set-up 
In all the experiments described here, use was 
made of a permanent magnet having a magnetic 
induction of 0.44 Wh/m?. The RF voltage was 
1 volt r.m.s. This voltage was taken from a Philips 
signal generator, type GM 2653. The electron current 
of the omegatron was stabilized. After amplifica- 
tion, the ion current was measured with a highly 
stable DC amplifier with negative feedback *). 
The smallest measurable current was 10-14 A. 
The resolution m/Am of the omegatron depends 
markedly on the magnitude of the magnetic induc- 
tion B, as appears from the equation: 
I Se TE) 
Am 2Epem 
3) J. Pelchowitch, Philips Res. Repts 9, 1, 1954. 
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Fig. 2. The omegatron with side plates in a glass envelope. 
The cathode is on the right, and on the extreme left is the 
collector for the electrons which have traversed the omegatron. 
The ions formed in the space inside the omegatron move to 
the ion collector. This is connected to the central metal pin 
fused into the glass tube and projecting below it. 


where S, is the distance from the point where the 
ion originated to the ion collector and Eye is the 


amplitude of the RF field. 


The resolution of the omegatron can easily be derived from 
the equation given in the literature for the radius r of the ion 
orbits *): 

Ene A 2 
r= sy, Sin of obit Semmre Gader oa C8) 
where ¢ is the difference between the resonance frequency ¢ 
of an ion having mass m and charge e and the frequency w 
of the applied RF field: ¢ = |w)— @,|. It is assumed that 
€ << We. The ions cannot pass beyond a circle of radius 


Send NUNIT 0 ge eee ete) 


Ifex< Ens/ BS), then rmax = Sp and ions of resonance frequency 
@®, reach the ion collector even though the frequency  ) of 
the RF field differs from we. If ¢ > Ent/ BS, Tmax <{ So and 
these ions do not arrive on the collector. It follows, then, that 
ions having a resonance frequency @¢ will continue to fall on 


*) H. Sommer, H. A. Thomas and J. A. Hipple, Phys. Rev. 
82, 697, 1951. 
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the collector while the frequency of the RF field is varied be- 


tween (, + € and My e. It will then only just be possible 
to discriminate completely between a mass m and a mass 
m-+- Am, provided the difference 4m, between the resonance 
frequencies of the two masses is equal to 2e. We can thus 


write the resolution as 


m Meo We 


eBS, 
2 Eng m ; 


6 
1m IWe 2e (6) 

In the experimental arrangement used the reso- 
lution is theoretically (i.e. according to equation (6)) 
such that two kinds of ions whose mass difference 


is Am = 1 can be completely separated if the 


mass number mis not greater than 34. In actual fact 
this was found to be the case up to m = 23. In order 
to improve the resolution, we have recently ob- 
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tained a permanent magnet whose magnetic induc- 
tion B is 0.92 Wh/m?. Fig. 3 shows a photograph of 
this magnet together with the other equipment used. 
The picture tube under investigation is held in an 
adjustable yoke. The distance between the poles 
of the magnet is 5.0 cm, and the diameter of the 
pole pieces is 13.0 em. At a distance of 2.0 em from 
the centre the variation in the magnetic induction is 
only 0.1%. Theoretically the resolution is now such 
that two kinds of ions with Am = 1 can be 
completely separated if m is not greater than 73. 
Experimentally, effective separation is found up 
to m = 49, which means that for example pro- 
pane gas (C;H,, mass number 44) can be fully 
distinguished from its “neighbours” using the large 


ie. 3. Experimental set-up, consisting of a permanent magnet (magnetic induction B — 
-92 Wh/m?), a DC amplifier for the ion current, a Philips signal generator, type GM 2653 
} 


which supplies the RF voltage, and ancillary 
the omegatron. The television picture tube ur 


equipment for the other voltages needed for 


ider investigation is connected to an omega- 


tron located between the pole pieces of the magnet. 
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magnet; but not when using the small magnet. 
Fig. 4 shows the spectrum of a mixture of propane 
and a trace of methane. This spectrum was recorded 
using the magnet with induction B = 0.44 Wb/m?. 


Fig. 4. Mass spectrogram of residual gases in a television pic- 
ture tube, showing ion current i+ of omegatron as a function 
of the preset mass number m. The tube contains propane 
C,H, and a small quantity of methane CH,. The RF voltage 
was 1 V. The magnet used for this recording has a magnetic 
induction B= 0.44 Wh/m?. 


Fig. 5 shows the same spectrum, but now with an 
induction B = 0.92 Wh/m?. The improved resolu- 
tion is particularly noticeable for the mass numbers 


43 and 44. 


Gas analyses on picture tubes in various stages of 


manufacture 


In the manufacture of television picture tubes 
it is important to know not only the composition 
of the residual gases in the finished tube, but also 
which parts of the tube contribute most to this 


£x10" 
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Fig. 5. The same mass spectrum as in fig. 4, but now using a 
magnet with B= 0.92 Wb/m*. The voltages applied to the 
omegatron were the same as in the previous analysis. The 
resolution is substantially improved. 
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residue, and what gases they contribute. This in- 
formation is obtained by studying unfinished tubes, 
i.e. tubes which have not yet been fitted with a 
fluorescent screen and/or other components. Ana- 
lyses are also carried out after the various operations 
following seal-off. 

In all cases an omegatron is fixed to the side of the 
tube immediately below the screen. First of all we 
investigate the glass bulb of a picture tube which 
contains only an unevaporated barium getter, but 
otherwise none of the “normal’’ components, such 
as fluorescent screen and electron gun. In addition 
to the omegatron a Philips ionization gauge, type 
53 EM °), is sealed to the tube. The only operation 
undergone by the tube is evacuation for two hours 
on an oil-diffusion pump. During this process 
the tube is heated for 20 minutes at 400 °C and 
allowed to cool to about 100 °C, after which 
it is sealed off. It is assumed that a more or 
less stationary state sets in 20 hours after seal- 
off. The gas composition is then measured with 
the omegatron. The gases found are hydrogen, me- 
thane, water vapour, carbon monoxide and carbon 
dioxide. Next, the getter is evaporated and a further 
measurement is made. The chemically active gases 
have now been absorbed by the getter, and in their 
place the tube contains a great deal of hydrocarbons, 
probably due to reactions between water vapour 
and carbon in the getter. The hydrocarbons are not 
taken up by the getter. A certain amount of argon, 
absorbed by the barium getter during its manufac- 
ture, is also found. The ionization gauge is now 
switched on for 15 minutes. The pumping action of 
this gauge depends on two effects. In the first place, 
hydrocarbons are decomposed at the hot cathode, 
which is at a temperature between 1500° and 2000°C. 
The fragments react with the getter, thereby reducing 
the pressure. Secondly, the electron current of 2.5 
mA produces ions in the gas, which are captured 
by the collector and thus removed from the gas. 
This too results in a pressure drop. At the end of the 
15 minutes the ion gauge is switched off and the 
tube is left for 14 days without anything being done 
to it. At the end of this period only a very slight 
increase in pressure is found. 

In fig. 6 the results of gas analyses done at four 
different times are compared: a) before and b) after 
evaporation of the getter, c) after 15 minutes of 
drawing current in the ionization gauge, and d) 
after 14 days. In the right-hand column the getter 
had been preheated (subjected to previous degassing) 
in the left-hand column it had not. The initial pres- 


5) E. Bouwmeester and N. Warmoltz, Philips tech. Rev. 
17, 121-125, 1955/56. 
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Fig. 6. Composition of the residual gases in unworked bulbs of 
television tubes. a) Before evaporation of barium getter; 
b) after evaporation of getter; c) after Philips 53 EM ionization 
gauge had drawn a current of 2.5 mA for 15 minutes; d) 14 
days later. In the left column the getter had not been preheated; 
in theright column it had. Seale for partial pressures in torrs 
(mm Hg), likewise in following figures. 


sure is lower when the getter has been preheated, 
but after 14 days no difference is to be found. 

Fig. 7 shows a comparison between four different 
bulbs analysed at four different phases in the pro- 
duction process, without the getter having been 
evaporated. The first bulb (a), as in the previous 
case, Contains no screen or components. The second 
(b) has been provided with a fluorescent screen, and 
in the third (c) the screen has been given the con- 
ductive coating of aluminium. The last bulb (d) 
also contains the electron gun. The cathode in the 
gun was degassed and activated when the bulb was 
still on the oil-diffusion pump. All four bulbs were 
subjected to the pumping process described above. 
Twenty hours after seal-off the gas composition was 
determined with the omegatron. It is noticeable 
that the hydrogen pressure is lower and the nitro- 
gen and argon pressures higher in the second bulb. 
The hydrogen has probably been taken up by the 
screen, and the higher nitrogen and argon pressure 
is presumably due to the uptake of these gases from 
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the atmosphere during the production of the screen. 
Bulb c¢ differs hardly at all from bulb b, but the 
finished tube with the electron gun has a lower 
nitrogen pressure and a higher carbon monoxide 
pressure. The carbon dioxide pressure has also gone 
up, probably as a result of the gas given off by the 
oxide cathode. 

After evaporation of the getter in the finished tube 
the gas composition changes radically, and in most 
tubes the residual gases consist solely of methane and 
argon. In some tubes, however, other hydrocarbons 
are found, but with a much lower partial pressure 
than methane. 

Fig. 8 shows how the composition of the residual 
gases in a finished picture tube changes during 
operation. Measurements were started twenty hours 


Sj SF S&S 
© N 
IQ 


ras 


OOO 


2 
"Was 


x 


Pas 


lo- 
2 


Z 
2S 


XO 


x 


K2 
bx 


S = Ss 
ra 
|o- 


Ar 


\N 4 

N Z 

ELE 
tO ° 
SPSLS 


rn 
PZ) 
CH, = 
N2 Rea 
SSS 


Ar 


5387 5388 


Fig. 7. Composition of residual gases in picture tubes in various 
stages of production. a) Unworked glass bulb; 6) bulb with 
fluorescent screen; c) bulb with fluorescent screen and conduc- 
tive aluminium coating; d) tube complete with electron gun. 
The getter was not evaporated in any of the tubes. 


Fig. 8. Composition of residual gases in a finished picture tube 
at four different times. a) 20 hours after seal-off, with the ca- 
thode cold; very little change occurs when the cathode is raised 
to working temperature. b) during the first hour of screen 
bombardment; c) after 250 hours of screen bombardment, 


measured during operation; d) immediately after switching 
the tube off. 
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after seal-off (a). While the cathode is heating up, 
very little change is found, but as soon as the screen 
is scanned by an electron beam, the nitrogen and 
hydrogen pressures increase and the hydrocarbons 
vanish (b). After scanning of the screen for 250 
hours with an electron current of 275 uA — the 
electron energy being 16 keV — the hydrogen and 
nitrogen pressures are found to have decreased (c). 
After scanning is stopped, the only residual gas 
which can be detected is argon (d). The pumping 
speed of the getter film is apparently not high 
enough to deal initially with the large production 
of hydrogen and nitrogen. We shall return later 
to the very sharp drop in the quantity of hy- 
drocarbons at the beginning of the scanning of 
the screen. 


Breaking open the picture tubes under vacuum 


We have seen how the composition of the 
residual gases in picture tubes can be determined 
by connecting an omegatron to the tube at the 
outset. It may happen, however, that a tube from 
normal production is not working properly because 
the gas pressure is too high. In such a case we want 
to be able to attach the omegatron to the rejected 
tube to ascertain the gas responsible for the fault. 

This is in the 


done 
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is made to break the thin glass wall, thus 
establishing communication between picture tube 
and omegatron. Fig. 9 shows a diagram of the 
whole arrangement. The part shown enclosed in the 
dotted lines is degassed for 15 hours at 400 °C 
before being connected to the rest of the system. 
At the same time the glass connecting tube is heated 
by hot wires to 200°C. During this process the 
epoxy resin does not get hotter than 100 °C, and 
the picture tube as a whole remains at room tem- 
perature. After 15 hours on the diffusion pump and 
outgassing, the pressure is reduced to a value be- 
tween 110-8 and 5x 10-° torr. The residual gases 
are then hydrogen, nitrogen, carbon dioxide and 
sometimes a small amount of water vapour. 

After the system has been cooled to room tem- 
perature, the high-vacuum valve K, is closed. This 
is a glass valve containing no grease and is closed by 
means of a magnet. Once the glass membrane be- 
tween the system and the picture tube has been 
broken, the gas content of the tube can be analysed 
with the omegatron. 


For this method to have any value the gas desorption of the 
epoxy resin must obviously be low. In the system without 
epoxy resin obtained by sealing off the glass tube above the 
connection to the picture tube, a pressure of 10~° torr can 


following way. First of all I 
an area 5 mm in diameter | O 

in the neck of the bulb near 0 
the electron gun is ground 
down with an ultrasonic 
drill to a thickness of only 
0.5 mm °). A thin glass | 
tube is cemented round the Pe 


weakened area with epoxy 


resin. The glass tube is 4 
connected to a pump 
system incorporating an 
omegatron, and a metal 
weight is introduced into 
the tube. The resin is hard- 
ened by heating at about R, 
100 °C for half an hour, 
resulting in a strong va- 
At a 
suitable the 
weight, which has a sharp 


steel point and can be 


cuum-tight joint. 
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Fig. 9. System for breaking open a television tube and analysing the residual gases. A 
glass tube is cemented with epoxy resin E round a spot in the neck of a picture tube T 
where the wall has been ground very thin. The connecting tube is degassed by a heater 
coil H. The part of the system inside the broken line can be outgassed in an oven at 400 °C. 
This part comprises an omegatron O, an ionization gauge I, a weight G and a glass high- 


vacuum valve K,. M, Penning gauge. K, second high-vacuum tap. K,-K, grease cocks. 


6) R. H. Collins and J. ©. 
Turnbull, Vacuum 10, 27, 
1960 (No. 1/2). 


M, Pirani (hot-wire) gauge. R, reservoir containing carbon monoxide. V expansion volume 
of 5.21. C liquid-nitrogen trap. L, capillary. The mercury diffusion pump is connected 
to P. The control valve K, and methane-filled bulb R, can be connected to point A. 
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easily be achieved. The residual gas is then primarily nitrogen. 
The lower pressure must therefore be attributed entirely to 
the absence of gas given off by the epoxy resin. Provided the 
pressure in the picture tube is much higher than 5x 10-° torr, 
which is always the case in practice, this method is therefore 
permissible. 

Reasonable agreement is found between the composition 
of the residual gases in the picture tube before it is broken open 
(measured with a second omegatron fixed directly to the bulb) 
and the composition found after the tube is broken open. 
The only difference found after breaking open the tube is that 
the pressures of the water vapour and carbon dioxide are lower, 
owing to absorption on the well-degassed components of the 


pump system. 


Fig. 10 shows the composition of the residual gases 
in two picture tubes that were rejected on the grounds 
of excessive gas pressure and were broken open in 
the manner described. These tubes were found to 
contain hydrocarbons, water vapour and argon. 
The hydrocarbons were probably due to reactions 
between water vapoar and carbon compounds in 
the getter film“). In other picture tubes the exces- 
sive pressure was sometimes found to be attribut- 
able to argon. The latter tubes showed no leakage, 
as was proved by storing them for a considerable 
period in plastic envelopes filled with helium. At 
the end of this time the omegatron registered no 
higher helium pressure. It is probable that air en- 
tered these tubes during the pumping process. The 
active gases would then be taken up by the getter 
film, but not the inert gas argon. 

The remainder of the system represented in fig. 9 
serves for investigating the gettering capacity of the 
barium getter in the picture tube 8). For this meas- 
urement K, is closed and K, opened. The carbon 
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Fig. 10. Composition of residual gases in two picture tubes from 
normal production, rejected because of their excessive gas 
pressure. 


‘) F. H. R. Almer, H. J. R. Perdijk and P. G. van Zanten, 
Advances in yacuum science and technology, Proc. Ist 
Int. congress on vacuum techniques, Namur, 10th-13th 

. aa ae Part 2, p. 676, Pergamon Press, Oxford 1960. 

*) J. J. B. Fransen and H. J. R. Perdijk, Philips tech. R fra 

19, 290, 1957/58. ; ie x 
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monoxide is admitted from a bulb having a 
known volume V of 5.21. The pressure p, in V is 
measured with a Pirani (hot-wire) gauge, which is 
specially calibrated for carbon monoxide. In this 
way a known quantity of gas can be admitted to 
the getter film. The pumping speed S of the getter 
can be found from the equation 


eee pee elit ese 


Pi 


GR) 


where p, is the carbon-monoxide pressure in the 
picture tube and L the conductance between the 
bulb and the picture tube. The conductance is almost 
entirely determined by the capillary Ly. The pres- 
sure p, in the picture tube can be measured by using 
the electron gun of the tube itself as an ionization 
The procedure is illustrated in fig. I1. 


gauge °). 
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Fig. 11. Pressure measurement in a picture tube using the elec- 
tron gun itself as an ionization gauge. K cathode at earth 
potential. G, grid negative to earth. G, positive grid which 
collects the electrons. The ion current, which is a measure of the 
pressure, is measured with G,. 


The cathode is kept at earth potential. With the 
aid of the first negative grid the electron beam is 
directed on to the second grid, which has a potential 
of 300 V positive to the cathode. All other electrodes 
are given a potential of 40 V negative to the cath- 
ode. These electrodes are used to determine the ion 
current, which in this case is a measure of the pres- 
sure. 

The total amount of carbon monoxide that a 
getter film in a picture tube can take up depends on 
the evaporated quantity of getter material. A re- 
presentative value is | torr-litre when 45 mg of 
barium is evaporated in a 21” picture tube. 

Other kinds of electron tubes that can be broken 
open and analysed in this way are oscilloscope tubes, 
transmitting tubes and television pick-up tubes, 
such as vidicons and image orthicons. Generally 


®) J. de Gier, Ingenieur 65, O 36, 1953. 
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speaking it can be said that any electron tube can 
be investigated in this way, provided the quantity 
of gas in the tube is greater than 10-7 torr l. 

An example of the analyses of an image orthicon 
can be seen in fig. 12. In this case the residual gases 
were analysed four months after manufacttre, and 
were found to consist of methane, ethane and argon 
from the barium getter. Also found was a partial 
pressure of 210-° torr for helium, which had 
diffused through the glass wall from the outside 
atmosphere. The quantity Q of helium that diffuses 
per unit time through a surface of area A and 
thickness d, at a pressure difference across the wall 


Ap, is given by: 
Q = ——_.. ..... (8) 


Here K is the permeability of helium of the boro- 
silicate glass used, given by '): 


torr-litre Xmm (thickness glass wall) 


K=1.5x10-% : =. 
em? (surface) x sec < torr (He pressure diff.) 


In the case of the image orthicon, A was 600 cm?, 
Ap was 3x10? torr, and the average wall thick- 
ness d 2.2 mm. The volume of the orthicon was 
740 cm, giving a calculated pressure rise after four 
months of 1.9 10~° torr. This is in excellent agree- 
ment with the measured helium pressure of 2 x 10~® 
torr. 


Decomposition of hydrocarbons in a picture tube 
during operation 

When the electron gun of a picture tube is used 
as an ionization gauge in the manner described 
above, the partial pressures of the hydrocarbons 
measured with the omegatron are found to decrease 
very rapidly. This is due to the disintegration of hy- 
drocarbons on the hot cathode, which is at a tempera- 
ture of 800 °C, and to the subsequent uptake of the 
fragments in the getter film. Furthermore, the colli- 
sion of electrons with gas molecules gives rise to ions 
which are captured by the negative electrodes. 
Both mechanisms contribute to the reduced partial 
pressures of the hydrocarbons. A plot of the de- 
crease measured in a picture tube is shown in fig. 13. 
The heavier hydrocarbon molecules decompose 
more rapidly than the light ones, so that the partial 
pressures of the heavier molecules decrease faster 
than that of methane molecules, for example, which 
are light. The pumping speed for the latter can be 
found from the formula: 


Vgc Spe Ry it Seen ) 


10) A, Klopfer and W. Ermrich, Vacuum 10, 128, 1960 (No. 1/2). 
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where V is the volume, p the pressure, dp/dt the 
pressure change per unit time and S the pumping 
speed. After integration we find: 
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Fig. 12. Composition of the resi- 
-8 dual gases in an image orthicon 
10 four months after manufac- 


ture. Note the high helium 
pressure. 
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where py is the pressure at t = 0. The pumping 
speed can thus be found from the slope of the 
curve obtained by plotting the logarithm of the 
pressure against time. The pumping speed found 
in this way for methane is roughly 0.004 l/sec. 

The pumping speed for hydrocarbons is much 
higher during normal operation of the tube, when a 
screen is bombarded by electrons of energy 16 keV. 
The exceptionally high pumping speed cannot be 
explained by the longer path of the electrons 
through the tube and the consequent larger number 
of ions formed. The production of ions does not even 
increase linearly with the length of the electron 
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Fig. 13. Decrease in the partial pressure of hydrocarbons in a 
picture tube. The logarithm of the pressure is plotted versus 
time in hours. Electron current 400 uA. The electron gun of 
the tube was used as an ionization gauge. 
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trajectories, since the electrons now have a higher 
energy and therefore their ionization probability 
is less than when the electron gun is used as 
an ionization gauge. The conclusion, then, is that 
the high pumping speed for hydrocarbons is 
attributable to the fluorescent screen. This is 
clearly demonstrated in fig. 14, where the loga- 
rithm of the methane pressure is plotted as a func- 
tion of time from the beginning of the screen bom- 
bardment. Curve | relates to a finished picture tube, 
and curve 2 to a tube without a fluorescent screen 
but containing the conductive aluminium coating. 
In the latter tube the pumping speed for methane 
is again 0.004 I/sec, compared with 2 I/sec in the 
finished tube. A similar comparison of two tubes, 
one of which had a silicate binder added to the 
fluorescent screen and the other not, revealed that 
the silicate binder is responsible for the high pump- 
ing speed for hydrocarbons. The form in which 
the silicate binder is present after preparation of 
the screen shows much resemblance to silica gel, 
a substance possessing a large internal surface and 
capable of adsorbing large amounts of gas. 

The fact that the fluorescent screen is responsible 
for the pumping effect is again clearly apparent 
from the effect of the accelerating voltage. The 
effect is not perceptible until 1.5 kV, i.e. the same 
potential at which the electrons are able to pene- 
trate through the aluminium layer to the phosphor 
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Fig. 14. Adsorption of methane during bombardment of the 
paecy in a finished picture tube (1) and in a tube without 
uores¢ ae 2 ee : 
uorescent screen (2). Constant pumping speeds are inferred 
from the slopes of the two curves. 
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coating. Up to 4 kV the pumping speed continues 
to rise, and above this value it is independent of 
the voltage. 

Various experiments were done to investigate the 
mechanism of this pumping effect. A finished pic- 
ture tube was connected by a glass tube to the pump- 
ing system shown in fig. 9. The connecting tube was 
temporarily closed off by a thin-walled glass bulb. 
The bulb R,, filled with methane, was connected 
to the pump system at 4 via a control valve K;. 
The part of the system shown within dotted lines 
in the figure was degassed for 15 hours at 400 °C. 
The picture tube was kept at room temperature. 
At the end of this period the total pressure measured 
with the omegatron was roughly 10-° torr. The 
residual gases consisted almost entirely of hydrogen, 
nitrogen and carbon dioxide. Communication with 
the picture tube was then established by breaking 
the thin-walled bulb with a steel ball. The high- 
vacuum valve K, was closed, after which the picture 
tube was given a methane pressure of about 5 x 10-6 
torr by adjustment of the control valve K,. The 
omegatron was adjusted so as to bring ions having 
a mass 16 into resonance, making it possible to 
determine the methane pressure continuously 
throughout the experiment. During this time the 
electron beam scanned a variable area on the fluores- 
cent screen. Fig. 15 shows the logarithm of the 
methane pressure plotted versus time from the 
moment of switching on the tube, for five different 
scanning areas. The current density of the electron 
beam at the screen was at all times 1 pA/em?, 
the accelerating voltage 16 kV. Whether the barium 
getter was evaporated or not had no effect on the 
results. 

Similar experimental runs were carried out for 
different current densities at the screen. They all 
gave results resembling those in fig. 15, except that 
the pressure decreased more rapidly the higher the 
current density. 

The characteristic shape of the curves in fig. 15 
may be explained qualitatively as follows. We sup- 
pose that the electrons, in bombarding the screen, 
give rise to active centres in the silicate binder of 
the fluorescent coating, as a result of which the 
incident methane molecules are trapped and at the 
same time hydrogen is desorbed. During the initial 
bombardment of the screen the number of centres 
increases and so does the pumping speed. Assuming 
that the number of centres reaches a maximum 
after a certain time, it follows that the pumping 
speed must then become constant. This is the case 
in the linear portion of the curve in fig. 15. After 
some considerable time a state of equilibrium is 
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Fig. 15. Variation of methane pressure in a picture tube. The 
logarithm of the pressure is plotted as a function of time in 
minutes for five different values of the scanned area. The 


current density of the electron beam was kept constant (j = 
1 wA/cm? on the screen). Va = 16 kV. 


reached in which the pressure remains constant. 
In this state the adsorption of methane is equal to 
its desorption. 


The mechanism may be described in somewhat greater detail 
as follows. Let ny be the maximum number of centres formed 
per unit area and NV the number of centres already produced 
on the scanned area A. The number of centres that can still 
be produced is then Anj—N. Experimeni shows that the rate 
of change of An,—N is proportional to the current density j 
at the screen. It is also found that the active centres have a 
certain lifetime and that the number vanishing per unit time 
is proportional to the number present. The number of centres 
that disappear as a result of methane adsorption is negligibly 
small. The number of centres appearing or vanishing per unit 
time is therefore given by: 


dN 


ore = k, j (An, — NN) — kN. eee (el) 
where k, and k, are constants. At equilibrium we have: 
dN 
— =0, 12 
dt (12) 
and therefore, putting k,/k, = k, we find: 
x ky jAng = WIS ; (13) 
kyj+k,  kyrl 


The number of methane molecules adsorbed by the screen per 
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unit time, i.e. V dp/dt, is proportional to the number of active 
centres and to the pressure: 


—- =—alNp, storey staal (Las) 


where «@ is a proportionality factor. Combining this expression 
with equation (9) it follows that the pumping speed S is equal 
to +alN, or, denoting the maximum pumping speed per unit 
area, ang, by so: 

kjAsy 
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In the linear portion of the “characteristic” in fig. 15 the 
pumping speed is constant and the number of centres formed 
under the given conditions is maximum. The maximum pump- 
ing speed per cm? is sy = 0.01 I/sec, and the constant k = 2.0. 
The formula derived shows good agreement with the experi- 
mental values of the pumping speed found for different values 
of scanned area A and current density j at the screen. The 
formula is valid up to a scanned area of 300 cm?, above which 
the increase in pumping speed is less than follows from the 
formula. 


The high pumping speed for methane in a picture 
tube during operation considerably reduces the 
partial pressure of the methane. A finished picture 
tube, when not in operation, often contains more 
methane than any other gas. This is due to the fact 
that methane may be produced at room tempera- 
ture by reactions between water vapour and carbon 
in the getter film. There is no uptake of methane 
in the barium getter itself. If the methane were not 
to disappear during operation, the ionic bombard- 
ment of the cathode would be very much stronger 
and the life of the tube would therefore be shortened. 
This illustrates clearly how important it is to know 
the composition of the residual gases and to inves- 
tigate the various processes that may take place 
inside a television picture tube. 


Summary. In the manufacture of television picture tubes it is 
important to know not only the total pressure in the tubes 
but also the composition of the residual gases. This information 
can readily be obtained with an omegatron. After briefly 
describing the principle of the omegatron for quantitative 
analyses and the experimental system used, the authors dis- 
cuss the desorption of gas by various parts of the picture tube. 
A method of analysing the composition of the gas in sealed-off 
tubes is described and some results are given. Finally, experi- 
ments are discussed which demonstrate the high uptake of 
methane by the fluorescent screen in a picture tube during 
electron bombardment. This removal of desorbed methane 
takes place continuously during the normal operation of a 
picture tube. It is attributable to the silicate binder used in 
the screen, and it guards the cathode against damage by ion 
bombardment. 
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2865: F. N. Hooge: Relation between electro- 
negativity and energy bandgap (Z. phys. 
Chemie Neue Folge (Frankfurt a.M.) 24, 
275-282, 1960, No. 3/4). 

It is shown that there exists a relation between 
the energy gap of crystals with formula AB and the 
electronegativities of the atoms A and B. The values 
of the electronegativities are the only quantities 
needed for the estimation of the width of the energy 
gap. 

2866: G. Thirup: The application of phase-locking 
techniques to the design of apparatus for 
measuring complex transfer functions (J. 
Brit. Instn. Radio Engrs. 20, 387-396, 1960, 
Noes): 

The theory of phase-lock synchronization is given 
briefly. Two devices for measuring complex voltage 
ratios are described. In the first, which has a fre- 
quency range 1-110 Mc/s, the voltage ratio to be 
measured is converted to a constant intermediate 
frequency of 450 ke/s. The frequency difference 
between the signal generator and the local oscillator 
is kept constant by coupling the tuning shafts 
mechanically and by a phase-lock synchronizing 
device; special means for extending the pull-in range, 
making use of a coarse and fine electronic adjust- 
ment, are described. The second instrument has a 
frequency range of 30-700 Mc/s, and here the local 
and signal generators are free-running and are 
separately tuned manually. A double frequency 
conversion system is used. The frequency of the 
second local oscillator, by means of a phase-lock 
device, is kept 450 ke/s below the frequency 
difference of the signal and first local oscillator. 
A description is also given of how the complex 
measuring results can be displayed on a cathode-ray 
tube. 


2867: M. Avinor: Gold-activated (Zn,Cd)S_phos- 
phors (J. Electrochem. Soc. 107, 608-611, 
1960, No. 7). 


Gold is shown to produce three emission bands in 


CdS, at 640, 800 and 1150 mu. The long-wave band 
appears only when a coactivator is used, while the 
short-wave bands are observed with activation by 
gold alone. The 1150 my band in CdS is shown to 
correspond to the 530 my gold band in ZnS. 


2868: H. Bremmer: The propagation over an in- 
homogeneous earth considered as a two- 
dimensional scattering problem (Electro- 
magnetic wave propagation, int. Conf. spon- 
sored by the Postal and Telecomm. Group 
of the Brussels Universal Exhibition, edited 
by M. Désirant and J. L. Michiels, pp. 
253-260, Academic Press, London 1960). 

This article deals with a two-dimensional integral 
equation for the distribution of field strength across 

a flat inhomogeneous earth. The equation in question 

is based on an approximative boundary condition, 

which accounts for the local distribution of the 
electrical constants. An infinite series representing 
the solution is interpreted in terms of scattering 
effects. The result applying to a special situation 

(two homogeneous regions separated by a straight 

boundary) is compared with that derived from the 

conventional one-dimensional equation which con- 
stitutes a saddlepoint approximation of the above 
two-dimensional equation. 


2869: C. Z. van Doorn and Y. Haven: Anisotrope 
kleurcentra in KCl (Ned. T. Natuurk. 26, 
216-220, 1960, No. 7). (Anisotropic colour 
centres in KCI; in Dutch.) 

When a KCl crystal is coloured by heating in K 
vapour (F absorption band at 5340 A) and then 
irradiated at room temperature with light absorbed 
in the F band, a new (M) band appears at 8080 A. 
After irradiation at 77 °K with polarized light ab- 
sorbed in the F band, both the F and M bands are 
found to be dichroic. The authors interpret this 
effect by assuming the M band to be caused by two 
neighbouring F centres. The merits of this model 
are compared with those of similar models proposed 


by Seitz and by Knox. 


